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Radioactive decay of nuclei via emission of a particles has been studied theoretically in the 
framework of a superasymmetric fission model using the double folding (DF) procedure for obtaining 
the a-nucleus interaction potential. The DF nuclear potential has been obtained by folding in the 
density distribution functions of the a nucleus and the daughter nucleus with a realistic effective 
interaction. The M3Y effective interaction has been used for calculating the nuclear interaction 
potential which has been supplemented by a zero-range pseudo-potential for exchange along with 
the density dependence. The nuclear microscopic a-nucleus potential thus obtained has been used 
along with the Coulomb interaction potential to calculate the action integral within the WKB 
approximation. This subsequently yields calculations for the half lives of a decays of nuclei. The 
density dependence and the exchange effects have not been found to be very significant. These 
calculations provide reasonable estimates for the lifetimes of a radioactivity of nuclei. 
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I. INTRODUCTION 



^ ' One of the first successes of quantum mechanics was the explanation of a-particle tunneling through barrier. About 
20 years after the first experimental observation of the a radioactivity [1,2], the theoretical explanation in terms of 
quantum mechanical barrier penetration [3,4] was provided which was, at least, qualitatively successful. Since then 
the a decay half life measurements with constantly improving techniques and the quest for theoretical explanation of 
their absolute magnitudes have continued unabated. The attractive features of a radiactivity are that the a particles 
| can be detected rather easily under favourable conditions such as high efficiency, low background and good energy 
i resolution. 

^5 About ten years ago, theoretical estimates of lifetimes for the a and various exotic decays have been provided in 

detail by the analytical superasymmetric fission model (ASAFM) [5,6] with reasonable success. This was followed 
1 1 ' by the cluster model (CM) calculations for the a decay half lives [7] for even more comprehensive database with 
similar success. But the nuclear interaction potentials used by both the theoretical approaches described above 
are phenomenological in nature and designed to fit the experimental data. The ASAFM uses a parabolic potential 
approximation for the nuclear interaction potential, which is a rather unusual a-nucleus interaction potential. The CM 
uses a strange cos-hyperbolic form for the nuclear interaction potential [7]. Nevertheless, they provide a benchmark 
against which more microscopically based treatments can be judged. 

In the present work, the nuclear potentials have been obtained microscopically by double folding the a and daughter 
nuclei density distributions with a realistic M3Y effective interaction. This is the ideal procedure of obtaining nuclear 
interaction energy for the a-nucleus interaction. Any liquid drop like properties such as surface energy are basically 
macroscopic manifestation of microscopic phenomena. A double folding potential obtained using M3Y effective 
interaction is more appropriate because of its microscopic nature. A potential energy surface is inherently embedded in 
this description. Moreover, the use of global microscopic nuclear potentials for a wide range of a-nucleus interaction 
is also theoretically a very thorough approach. The double folding potential has then been utilised within the 
superasymmetric fission model description. This superasymmetric fission model using microscopic potentials has 
been used to provide estimates of a decay half lives for a wide range of a-emitters. The semirealistic explicit density 
dependence on the M3Y effective interaction has also been employed to study the effects of density dependence on the 
a decay lifetimes. The penetrability of the pre-scission part of the potential barrier provides the a cluster preformation 
probability within a superasymmetric fission model [6] . 
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II. FORMALISM 



The microscopic nuclear potentials Vn(R) have been obtained by double folding in the densities of the fragments 
a and daughter nuclei with the finite range realistic M3Y effective interacion as 

V N (R) = J J P i(ri)p 2 (ri)v[\ri - f x + R\]d 3 ri d 3 r 2 (1) 

where the density distribution function p\ for the a particle has the Gaussian form 

pi{r) = 0.4229ea;p(-0.7024r 2 ) (2) 

whose volume integral is equal to A a {= 4), the mass number of a-particle. The density distribution function p 2 used 
for the residual cluster, the daughter nucleus, has been chosen to be of the spherically symmetric form given by 

p 2 (r)=po/[l + exp((r-c)/a)] (3) 

where 

c = r p (l-7c 2 a 2 /3r 2 p ), r p = 1.13A 1 / 3 and a = 0.54 fm (4) 

and the value of po is fixed by equating the volume integral of the density distribution function to the mass number 
Ad of the residual daughter nucleus. The distance s between any two nucleons belonging to the residual daughter 
nucleus and the emitted a nucleus is given by 

s=\f* 2 -f{+R\ (5) 

while the interaction potential between any such two nucleons v(s) appearing in eqn.(l) is given by the M3Y effective 
interaction. The total interaction energy E{R) between the a nucleus and the residual daughter nucleus is equal to 
the sum of the nuclear interaction energy, the Coulomb interaction energy and the centrifugal barrier. Thus 

E(R) = V N (R) + Vc(R) + h 2 l{l + l)/(2pR 2 ) (6) 

where fi = mA a A d /A is the reduced mass, A is the mass number of the parent nucleus and m is the nucleon mass 
measured in the units of MeV/c 2 . Assuming spherical charge distribution (SCD) for the residual daughter nucleus 
and considering the a-paticle to be a point charge, the a-nucleus Coulomb interaction potential Vc(R) is given by 

V C (R) = Z a Z d e 2 /R for R > R c 

= (Z a Z d e 2 /2R c ).[3 - (R/R c ) 2 } for R < R c (7) 

where Z a and Z d are the atomic numbers of the a-particle and the daughter nucleus respectively. The touching radial 
separation R c between the a-particle and the daughter nucleus is given by R c = c a + c d where c a and c d has been 
obtained using eqn.(4). If point charge distribution (PCD) is also assumed for the residual daughter nucleus then the 
a-nucleus Coulomb interaction potential has the simple form of Vc(R) — Z a Z d e 2 /R for all R. The energetics allow 
spontaneous emission of a-particles only if the released energy 

Q = M-{M a + M d ) (8) 

is a positive quantity, where M, M a and M d are the atomic masses of the parent nucleus, the emitted a-particle and 
the residual daughter nucleus, respectively, expressed in the units of energy. It is important to mention here that the 
correctness of predictions for possible decay modes, therefore, rests on the accuracy of the ground state masses of 
nuclei. 

In the present work, the half life of the parent nucleus against the split into an a and a daughter is calculated using 
the WKB barrier penetration probability. The assault frequecy v is obtained from the zero point vibration energy 
E v = {1/2)Tiuli = {l/2)hv. The half life T of the parent nucleus {A,Z) against its split into an a (A a ,Z a ) and a 
daughter (A d , Z d ) is given by 

T=[(/iln2)/(2^)][l + exp(^)] (9) 
where the action integral K within the WKB approximation is given by 
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K=(2/h) [2p(E(R)-E v -Q)] 1/2 dR (10) 

JRa 

where R a and Rb are the two turning points of the WKB action integral determined from the equations 

E(R a ) =Q + E V = E(R h ) (11) 

For the present calculations, from a fit to a selected set of experimental data on a emitters the following law [8] was 
found for the zero point vibration energies for the a decays 

£^ = 0.1045. Q for even(Z) — even(N) parent nuclei 

= 0.0962.Q for odd(Z) — even(N) parent nuclei 

= 0.0907.Q for even(Z) — odd(N) parent nuclei 

= 0.0767.Q for odd(Z) - odd(N) parent nuclei (12) 

which includes the shell and pairing effects. The values of the proportionality constants of E v with Q is the largest for 
even-even parent and the smallest for the odd-odd one. If all other conditions are the same one may observe that with 
greater value of E v , the life time is shortened indicating higher emission rate. The shell effects of a radioactivity is 
implicitly contained in the zero point vibration energy due to its proportionality with the Q value, which is maximum 
when the daughter nucleus has a magic number of neutrons and protons. 



III. CALCULATIONS 



The two turning points of the action integral given by eqn.(10) have been obtained by solving eqns.(ll) using the 
microscopic double folding potential given by eqn.(l) along with the Coulomb potential given by eqn.(7) and the 
centrifugal barrier. Then the WKB action integral between these two turning points has been evaluated numerically 
using eqn.(l), eqn.(6), cqn.(7), eqn.(8) and eqn.(12). The calculations have been performed using v(s), inside the 
integral of eqn.(l) for the double folding (DF) procedure, as only the M3Y effective [9] interaction 

v(s) = 7999exp(-4.s)/(4.s) - 2134exp(-2.5s)/(2.5s) (13) 

This interaction is based upon a realistic G-matrix. Since the G-matrix was constructed in an oscillator representation, 
it is effectively an average over a range of nuclear densities and therefore the M3Y has no explicit density dependence. 
For the same reason there is also an average over energy and the M3Y has no explicit energy dependence either. The 
only energy dependent effects that arises from its use is a rather weak one contained in an approximate treatment of 
single- nucleon knock-on exchange. The success of the extensive analysis [10,11] indicates that these two averages are 
adequate for the real part of the optical potential for heavy ions at energies per nucleon of < 20MeV. However, it is 
important to consider the density and energy dependence explicitly for the analysis of a-particle scattering at higher 
energies (> lOOMe^) where the effects of a nuclear rainbow are seen and hence the scattering becomes sensitive to 
the potential at small radii. Such cases were studied introducing suitable and semirealistic explicit density dependence 
[12,13] into the M3Y interaction which was then called the DDM3Y and was very successful for interpreting consistently 
the high energy elastic a scattering data. 

The entire calculations have been redone again using the density dependent M3Y effective interaction (DDM3Y) 
supplemented by a zero-range pseudo potential, representing the single nucleon exchange term [13]. In DDM3Y 
the effective nucleon-nucleon interaction v(s) is assumed to be density and energy dependent and therefore becomes 
functions of density and energy and is generally written as 

v(s, Pl ,p 2 , E) = t M3Y (s, E)g( Pl ,p 2 , E) (14) 

where t M3Y is the same M3Y interaction given by eqn.(13) but supplemented by a zero range pseudo-potential 

t M3Y = 7999exp(-4.s)/(4.s) - 2134cxp(-2.5s)/(2.5s) + J 0Q (E)5(s) (15) 

where the zero-range pseudo-potential representing the single-nucleon exchange term is given by 

Joo(-E) = -276(1 - 0.005E/A a ){MeV.fm 3 ) (16) 

and the density dependent part has been taken to be [13] 
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g( Pl ,p 2 , E) = C(l - /?(£)p? /3 )(l - (3(E)p 2 2 /3 ) 



(17) 



which takes care of the higher order exchange effects and the Pauli blocking effects. 

The energy E appearing in the above equations is the energy measured in the centre of mass of the a - daughter 
nucleus system and for the a decay process it is equal to the released energy Q. Since the released energies involved 
in the a decay processes are very small compared to the energies involved in high energy a scattering, (3(E) has 
been considered as a constant and independent of energy. The zero-range pseudo-potential Joo (E) is also practically 
independent of energy for the a decay processes and can be taken as — 276MeV.fm 3 . It is important to mention 
here that use of the nuclear interaction potential for exotic cluster decays obtained by double folding an effective 
nucleon-nucleon interaction with their respective densities was suggested in reference [14]. After the WKB action 
integral K, given by eqn.(10) has been evaluated, the half lives of the a decays have been calculated using eqn.(9) 
and eqn.(12). 

IV. RESULTS AND DISCUSSIONS 

For the present illustrative calculations, the same set of experimental data of reference [15] along with the rest from 
reference [16] for the a decay half lives have been chosen for comparison with the present theoretical calculations 
for even-even parent nuclei for which the experimental ground state masses for the parent and daughter nuclei are 
available. This set was selected because there is no uncertainty in the determination of the released energy Q (given 
by eqn.(8)) which is one of the crucial quantity for quantitative prediction of decay half lives and all the the parent and 
daughter nuclei have zero spins and positive parities. For the calculations using DDM3Y with zero range exchange 
also, the value of the normalization C appearing in eqn.(17) has been fixed at 1.0. The density dependent parameter 
13(E), which is supposed to be dependent on energy, has been kept constant and independent of Q. Optimum fit to 
the data has been obtained using /3(E)— 1.6. 

In Fig. 1 the experimental data for logarithmic a decay half lives [15,16] have been plotted against the mass numbers 
of parent nuclei along with the results of the present calculations for zero angular momentum of the fragments. In 
the figure the open circles depict the experimental data while the continuous line with solid circles represents the 
corresponding calculations of the present model using the M3Y effective interaction and the dotted line represents the 
present calculations using DDM3Y effective interaction supplemented by a zero-range pseudo-potential. The decay 
modes and the experimental values for their half lives have been presented in Table 1. The corresponding results of 
the present calculations of superasymmetric fission model with microscopic potentials are also presented along with 
the results of the liquid drop model (LDM) calculations of 2001 so as to facilitate the comparison of the results of the 
phenomcnological calculations with the present one. The exact Q values calculated by eqn.(8) using the experimental 
ground state masses [17] and used by the present calculations have also been presented. 

The results of the present calculations with M3Y or DDM3Y with pseudopotential have been found to predict the 
general trend of experimental data very well. The quantitative agreement with experimental data is excellent. The 
degree of reliability of the present estimates for the a decay lifetimes are better than the liquid drop description 
[15] for most of the cases. The recent result of the preformed cluster model [18] calculation for the a decay half-life 
of 242 Cm is also much worse compared to that obtained from present calculations. The degree of reliability of the 
present estimates for the a decay lifetimes is better than the recent results of the LDM calculations [19] presented in 
the table below. 



TABLE I. Comparison between Calculated a-decay Half-Lives using Point (PCD) and Spherical charge distributions (SCD) 
respectively for the Coulomb interaction and using effective interactions of M3Y and DDM3Y with zero-range pseudo-potential 
for the nuclear interaction. 
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logioT(s) 


Q(MeV) 


logioT(s) 


logioT(s) 


logwT(s) 


88 


222 


1.72 


6.68 


1.39(1.50) 


1.43(1.52) 


1.58 


88 


224 


5.74 


5.79 


5.34(5.41) 


5.39(5.44) 


5.50 


88 


226 


10.98 


4.87 


10.56(10.60) 


10.62(10.63) 


10.703 


90 


228 


8.07 


5.53 


7.71(7.75) 


7.77(7.78) 


7.781 


90 


230 


12.75 


4.78 


12.26(12.28) 


12.34(12.31) 


12.376 


92 


232 


9.69 


5.42 


9.34(9.36) 


9.42(9.39) 


9.337 


92 


234 


13.21 


4.86 


12.84(12.83) 


12.91(12.87) 


12.889 


94 


236 


8.01 


5.87 


7.87(7.87) 


7.96(7.91) 


7.954 
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94 


238 


9.54 


5.60 


96 


242 


7.11 


6.22 


90 


226 


3.49 


6.46 


90 


232 


18.18 


4.08 


92 


230 


6.54 


6.00 


92 


236 


15.23 


4.58 


94 


240 


13.39 


5.26 



(a) Taken from reference [16]. 



9.30(9.31) 


9.39(9.35) 


9.4423 


7.01(7.02) 


7.11(7.06) 


7.1485 


3.19(3.27) 


3.24(3.30) 


3.39 a 


17.63(17.62) 


17.72(17.66) 


17.76 a 


6.26(6.30) 


6.32(6.33) 


6.43 a 


14.80(14.78) 


14.89(14.82) 


14.99 a 


11.28(11.27) 


11.38(11.30) 


11.45 a 
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Above results show that the differences in the a-decay lifetimes obtained using the spherical charge distributions or 
point charge distributions for calculating the Coulomb interaction are small. Introducing the density-dependence and 
the zero-range pseudopotential do not alter the general trend of the results and also do not improve the quality of fit 
to the experimental data significantly. Use of the M3Y effective interaction alone for calculating the double folding 
nuclear interaction potentials is sufficient to provide reasonable estimates for the a-decay lifetimes. 



V. CONCLUSIONS 



The half lives for a-radioactivity have been analyzed with microscopic nuclear potentials obtained by the double 
folding pocedurc using M3Y and DDM3Y effective interactions respectively. This procedure of obtaining nuclear 
interaction potentials is based on profound theoretical basis. The results of the present calculations using M3Y or 
DDM3Y supplemented by a pseudo-potential arc in good agreement over a wide range of experimental data. Present 
calculations show that the differences in the results of the a decay lifetimes obtained using the spherical charge 
distributions or point charge distributions for calculating the Coulomb interaction are small. Use of the density- 
dependence and the zero-range pseudopotential also neither alters the general trend of the results nor improves 
significantly the quality of fit to the experimental data. Refinements such as introduction of dissipation while tunneling 
through the barrier or incorporating the dynamic shape deformations in the density distributions of the clusters may 
further improve results. It is worthwhile to mention that using the realistic microscopic nuclear interaction potentials, 
the results obtained for the a radioactive decay lifetimes are noteworthy and are comparable to the best available 
theoretical calculations. Such calculations may be extended to provide reasonable estimates of the lifetimes of nuclear 
decays by a emissions for the entire domain of exotic nuclei. 
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FIG. 1. Logarithmic a-radioactivity half-lives plotted against the parent mass number. The experimental data have been 
shown by the open circles. The solid circles represent the results of present calculations, using M3Y effective interaction, 
for the corresponding experimental data . The continuous line connects these calculated values. The dotted line (almost 
indistinguishable from the continuous line) represent the same calculations using DDM3Y with a zero-range pseudo-potential 
for the nucleon-nucleon effective interaction. For both the theoretical calculations spherical charge distributions for the Coulmb 
interaction has been assumed. 
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